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Abstract 

In a 1990 paper Helton and Young showed that under certain condi- 
tions the optimal solution of the Nehari problem corresponding to a finite 
rank Hankel operator with scalar entries can be efficiently approximated 
by certain functions defined in terms of finite dimensional restrictions of 
the Hankel operator. In this paper it is shown that these approximants 
appear as optimal solutions to restricted Nehari problems. The latter 
problems can be solved using relaxed commutant lifting theory. This ob- 
servation is used to extent the Helton and Young approximation result to a 
matrix-valued setting. As in the Helton and Young paper the rate of con- 
vergence depends on the choice of the initial space in the approximation 
scheme. 

1 Introduction 

Since the 1980s, the Nehari problem played an important role in system and con- 
trol theory, in particular, in the 7?°°-control solutions to sensitivity minimiza- 
tion and robust stabilization, cf., [9... In system and control theory the Nehari 
problem appears mostly as a distance problem: Given G in L°°, determine the 
distance of G to H°°, that is, find the quantity d := inf{||G - | F G H°°} 
and, if possible, find an F G H°° for which this infimum is attained. Here all 
functions are complex-valued functions on the unit circle T. It is well-known 
that the solution to this problem is determined by the Hankel operator H which 
maps H 2 into K 2 = L 2 Q H 2 according to the rule Hf = P-(Gf), where P_ is 
the orthogonal projection of L 2 onto K 2 . Note that H is uniquely determined 
by the Fourier coefficients of G with negative index. Its operator norm deter- 
mines the minimal distance. In fact, d — \\H\\ and the infimum is attained. 
Furthermore, if H has a maximizing vector ip, that is, <p is a non-zero function 
in H 2 such that = \\H\\ \\<p\\, then the AAK theory pQ [3] (see also [IS] ) 

tells us that the best approximation G of G in H°° is unique and is given by 

8(e«) = G(e«) - »p a.e. (1.1) 

By now the connection between the Nehari problem and Hankel operators is 
well established, also for matrix-valued and operator-valued functions, and has 
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been put into the larger setting of metric constrained interpolation problems, 
see, for example, the books (6] Chapter IX], [13l Chapter XXXV], Chapter 
I], [T7J Chapter 5] and [3J Chapter 7], and the references therein. 

The present paper is inspired by Helton- Young [TJ] . Note that formula (|l.ip 
and the maximizing vector <p, may be hard to compute, especially if H has 
large or infinite rank. Therefore, to approximate the optimal solution 
Helton- Young [14] replaces H by the restriction H = H\^Q Z n qH 2 to arrive at 

nt it\ r<i it\ {Hip)(e lt ) . . 

G(e ) = G(e ) — , a.e. (1.2) 

as an approximant of G. Here n is a positive integer, q is a polynomial and (p is 
a maximizing vector of H. Note that a maximizing vector of H always exists, 
since ranki? < n + degg, irrespectively of the rank of being finite, or not. 

In |14] it is shown that G is a computationally efficient approximation of 
the optimal solution G when the zeros of the polynomial q are close to the 
poles of G in the open unit disk D that are close to the unit circle T. To be 
more precise, it is shown that if G is rational, i.e., ranki? < oo, and \\H\\ is 
a simple singular value of H, then \\G — G\\oo converges to as n — > oo. This 
convergence is proportional to r n if the poles of G in ED are within the disc 
ES> r = {zGC||z|<r}, and the rate of convergence can be improved by an 
appropriate choice of the polynomial q. 

It is well-known that the Nehari problem fits in the commutant lifting frame- 
work, and that the solution formula (|l.ll) follows as a corollary of the commutant 
lifting theorem. We shall see that the same holds true for formula (|1.2[) provided 
one uses the relaxed commutant lifting framework of 8 ; cf., Corollary 2.5 in 

To make the connection with relaxed commutant lifting more precise, define 
R n to be the orthogonal projection of H 2 onto H 2 z n ~ 1 qH 2 , and put Q n = 
SR n , where S is the forward shift on H 2 . Then the operators R n and Q n both 
map H 2 into H 2 z n qH 2 , and the restriction operator H n := H\jf2 eznqH 2 
satisfies the intertwining relation V-H n R n = H n Q n . Here V- is the compression 
of the forward shift V on L 2 to K 2 . Given this intertwining relation, the relaxed 
commutant lifting theorem ,8, Theorem 1.1] tells us that there exists an operator 
B n from H 2 z n qH 2 into L 2 such that 

P-B n = H n , VB n R n = BQ n , \\B n \\ = \\H n \\. (1.3) 

The second identity in (ll.3[) implies (see Lemma \2 . 2 1 below) that for a solution 
B n to (jl.3|) there exists a unique function $„ e L 2 such that the action of B n 
is given by 

(B n h)(e lt ) = <$> n {e lt )h(e u ) a.e. (h e H 2 z n qH 2 ). (1.4) 

Furthermore, since H n has finite rank, there exists only one solution B n to (|1.3p 
(see Proposition 12.31 below) , and if ip n = tp is a maximizing vector of H n , then 
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this unique solution is given by (|1.4[) with $„ equal to 

it (g n V> n )(e ri ) (Hy)(e jt ) n 
$ " (e } = ^(e«) = ~m~' (L5) 

Thus G — G appears as an optimal solution to a relaxed commutant lifting 
problem. 

This observation together with the relaxed commutant lifting theory devel- 
oped in the last decade, enabled us to extent the Helton- Young convergence 
result for optimal solutions in [14] to a matrix-valued setting, that is, to derive 
an analogous convergence result for optimal solutions to matrix- valued Nehari 
problems; see Theorem 13.11 below . A complication in this endeavor is that for- 
mula (jl.lj) generalizes to the vector- valued case, but not to the matrix- valued 
case. Furthermore, in the matrix-valued case there is in general no unique so- 
lution. We overcome the latter complication by only considering the central 
solutions which satisfy an additional maximum entropy-like condition. On the 
way we also derive explicit state space formulas for optimal solutions to the 
classical and restricted Nehari problem assuming that the Hankel operator is 
of finite rank and satisfies an appropriate condition on the space spanned by 
its maximizing vectors. These state space formulas play an essential role in the 
proof of the convergence theorem. 

This paper consists of 6 sections including the present introduction. In 
Section [2] which has a preliminary character, we introduce a restricted version 
of the matrix- valued Nehari problem, and use relaxed commutant lifting theory 
to show that it always has an optimal solution. Furthermore, again using relaxed 
commutant lifting theory, we derive a formula for the (unique) central optimal 
solution. In Section [3] we state our main convergence result. In Section @] 
the formula for the (unique) central optimal solution derived in Section [2] is 
developed further, and in Section [5] this formula is specified for the classical 
Nehari problem. Using these formulas Section [6] presents the proof of the main 
convergence theorem. 

Notation and terminology. We conclude this introduction with a few words 
about notation and terminology. Given p, q in N, the set of positive integers, we 
write L 2 xp for the space of all gxp-matrices with entries in L 2 , the Lebesgue 
space of square integrable functions on the unit circle. Analogously, we write 
H 2 xp for the space of all gxp-matrices with entries in the classical Hardy space 
H , and K 2 X stands for the space of all gxp-matrices with entries in the 
space K 2 — L Q H 2 , the orthogonal compliment of H 2 in L 2 . Note that each 
F E L 2 xp can be written uniquely as a sum F = F + + F- with F + £ H 2 xp 
and F_ 6 K 2 xp . We shall refer to F + as the analytic part of F and to F_ 
as its co- analytic part. When there is only one column we simply write L 2 , 
H 2 and K 2 instead of L pxl , H pxl and K pxl . Note that L 2 , H 2 and K 2 are 
Hilbert spaces and K 2 = L 2 QH 2 . Finally, L^ xp stands for the space of all qxp- 
matrices whose entries are essentially bounded on the unit circle with respect to 
the Lebesque measure, and H^ p stands for the space of all gxp-matrices whose 
entries are analytic and uniformly bounded on the open unit disc U. Note that 
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each F £ L^ xp belongs to L 2 xp and hence the analytic part F + and the co- 
analytic part F_ of F are well defined. These functions belong to L 2 xp and 
it may happen that neither F + nor F_ belong to L"^ xp . In the sequel we shall 
need the following embedding and projection operators: 

E : C p -> H%, Eu{\) =u (z £ B); (1.6) 

U:K 2 q ^a, nf = ±-JJe- it f(e it )dt. (1.7) 

Throughout G G and ff : JZ? — > K 2 is the Hankel operator defined by 

the co-analytic part of G, that is, Hf = P_(G/) for each / G H p . Here P_ is 
the orthogonal projection of L 2 onto K 2 . Note that V-H = HS, where S is the 
forward shift on H p and V- is the compression to K 2 of the forward shift V on 

Finally, we associate with the Hankel operator H two auxiliary operators 
involving the closure of its range, i.e., the space X = ImH, as follows: 

Z-.X^X, Z = V-\ X , (1-8) 
W:H 2 ^X, Wf = Hf (feH 2 ). (1.9) 

Note that X := ImH is a ^--invariant subspace of K 2 . Hence Z is a well- 
defined contraction. Furthermore, if rankff is finite, then the spectral radius 
r spec (Z) is strictly less than one and the co-analytic part G_ of G is the rational 
matrix function given by 

G_(A) = (n|*)(AJ - Z)~ X WE. 

In system theory the right hand side of the above identity is known as the 
restricted backward shift realization of G_; see, for example, Section 7.1]. 
This realization is minimal, and hence the eigenvalues of Z coincide with the 
poles of G_ in D. In particular, r spcc (Z) < 1. Since V^H = HS, we have 
ZW = WS. Furthermore, Kcr H* = K 2 9 X. 



2 Restricted Nehari problems and relaxed corn- 
mutant lifting 

In this section we introduce a restricted version of the Nehari problem, and 
we prove that it is equivalent to a certain relaxed commutant lifting problem. 
Throughout M. is a subspace of H 2 such that 

S*MCM, Ker S* G M. (2.1) 

With M we associate operators Rm and Qm acting on H 2 , both mapping H 2 
into Ai. By definition Rj^ is the orthogonal projection of H 2 onto S*ftA and 
Qm — SRm- 
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We begin by introducing the notion of an .M-norm. We say that $ £ L qxp 
has a finite M-norm if S L^ for each h £ M. and the map h M> $/i is a 
bounded linear operator, and in that case we define 

\\$\\m =sup{||*ft|| z j \heM, \\h\\ H 2 < 1}. 

If M. is finite dimensional, then each $ £ L 2 „ has a finite .M-norm. Further- 
more, $ £ L^ xp has a finite .M-norm for every choice of A4, and in this case 
H^Ha-i < ll^lloo, with equality if M = H%. Note that $ £ L qxp has a finite 
.M-norm and G £ L^ xp imply G — $ has a finite .M-norm. 

We are now ready to formulate the .M-restricted Nehari problem. Given 
G £ L^ xp and a subspace M of H p , we define the optimal AA -restricted Nehari 
problem to be the problem of determining the quantity 

d M := inf{||G - F\\ M | F £ H qxp and F has a finite .M-norm}, (2.2) 

and, if possible, to find a function F 6 H qxp of finite .M-norm at which the 
infimum is attained. In this function F attaining the infimum is called 

an optimal solution. The suboptimal variant of the problem allows the norm 
F\\m to be larger than the infimum. When M = H p , the problem coincides 
with the classical matrix- valued Nehari problem in L^ p . In [15j [16] the case 
where M = H p S k H p , with k £ N, was considered. 

Proposition 2.1. Let G £ L^ xp , and let M. be a subspace of H p satisfying 
the conditions in (I2.1[) . Then the A4-restricted Nehari problem has an optimal 
solution and the quantity dj^t in (12.21) is equal to jm := where H : 

H p — > K q is the Hankel operator defined by the co-analytic part of G. 

We shall derive the above result as a corollary to the relaxed commutant 
lifting theorem [8] Theorem 1.1], in a way similar to the way one proves the 
Nehari theorem using the classical commutant lifting theorem (see, for example, 
[HI Section II.3]). For this purpose we need the following notion. We say that 
an operator B from M. into L q is defined by a $ £ L qxp if the action of B is 
given by 

(Bh)(e u ) = <5>{e lt )h{e lt ) a.e. {h£M). (2.3) 

In that case, $ has a finite .M-norm, and || || ^vi = When (|2.3|) holds we 

refer to $ as the defining function of B. The following lemma characterizes 
operators B from M. into L q defined by a function $ £ L qxp in terms of an 
intertwining relation. 

Lemma 2.2. Let M. be a subspace of H p satisfying (|2.ip . and let B be a bounded 
operator from Ai into L q . Then B is defined by a $ £ L qxp if and only if 
B satisfies the intertwining relation VBRj^ = BQm- I n that case, <fr(-)u = 
BEu(-) for any u £ C p and \\B\\ = \\®\\ M 

Proof. This result follows by a modification of the proof of Lemma 3.2 in [TTj . 
We omit the details. □ 
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Proof of Proposition I2TT1 Put = ||-^|m||- Recall that the Hankel opera- 
tor H satihes the intertwining relation V-H = HS. This implies V~H\mRm = 
H\mQm- Here Rm and Qm ar e the operators defined in the first paragraph 
of the present section. Since Q* m Qm — R%[Rm an d V is an isometric lifting 
of V- , the quintet 

{H\ m ,V-,V,R m ,Qm,Jm} (2.4) 

is a lifting data set in the sense of Section 1 in [5]. Thus Theorem 1.1 in [8 
guarantees the existence of an operator B from A4 into L 2 , with the properties 

P-B = H\ M , VBR M =BQ M , \\B\\= 7M . (2.5) 

By Lemma l2~2l the second equality in (|2.5[) tells us there exists a $ G L 2 xp 
defining B, that is, the action of B is given by (I2.3[) . As $(-)u = BEu(-), the 
first identity in (12. 5j) shows that G_ = and hence F := G - $ G H^ xp . 
Furthermore, 

||G-F|| A< = ||$|| A< = 11^11=7^, 

because of the third identity in (12.51) . Thus the quantity d,M in (12.21) is less than 
or equal to jm- 

It remains to prove that (1m > 1m- In order to do this, let F G H^ xp and 
have a finite A4-norm. Put $ = G — F. Then $ has a finite A^-norm. Let B be 
the operator from A4 into defined by $. Since i 7, G H^ xp , we have G_ = 
and hence the first identity in (|2.5[) holds with i? in place of i3. It follows that 

\\G-F\\ M = \\$\\ M = \\B\\ > \\H\ M \\ =1M- 

This completes the proof. □ 



In the scalar case, or more generally in the case when p = 1, the optimal 
solution is unique. Moreover this unique solution is given by a formula analogous 
to (|1.2p ; cfi, pp. This is the contents of the next proposition which is proved 
in much the same way as the corresponding result for the Nehari problem. We 
omit the details. 

Proposition 2.3. Assume p = 1, that is, G G and A4 a subspace of H 2 
satisfying (|2.1[) . Assume that H\m has a maximizing vector ip G M. Then 
there exists only one optimal solution F to the M.-restricted Nehari problem 
(|2.5[) . and this solution is given by 

F(e«) = G{e*) - a.e. (2.6) 

In general, if p > 1 the optimal solution is not unique. To deal with this 
non- uniqueness, we shall single out a particular optimal solution. 

First note that the proof of Proposition 12 . 1 1 shows that there is a one-to-one 
correspondence between the optimal solutions of the A^-restricted Nehari prob- 
lem of G and all interpolants for H\m with respect to the lifting data set (|2.4[) . 
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that is, all operators B from j\4 into L 2 satisfying (|2.5p . This correspondence 
is given by 



B H» F = G — $, where $ is the defining function of B. 



(2.7) 



Next we use that the relaxed commutant lifting theory tells us that among 
all interpolants for H\m with respect to the lifting data set (|2.4p there is a 
particular one, which is called the central interpolant for H\m with respect to 
the lifting data set (|2.4[) ; see Section 4]. This central interpolant is uniquely 
determined by a maximum entropy principle (see [5J Section 8]) and given by 
an explicit formula using the operators appearing in the lifting data set. 

Using the correspondence (|2.7p we say that an optimal solution F of the 
.M-restricted Nehari problem of G is the central optimal solution whenever 
$ := G — F is the defining function of the central interpolant B for H\m with 
respect to the lifting data set (12.4[) . Furthermore, using the formula given in 
[HI Section 4] for the central interpolant the correspondence (|2.7p allows us to 
derive a formula for the central optimal solution. To state this formula we need 
to make some preparations. 

As before j M = \\H\m\\- Not e that \\HP M S\\ < \\HP M \\ = \\H\ M \\, where 
Pm is the orthogonal projection of H 2 (C P ) on Ai. This allows us to define the 
following defect operators acting on H 2 (C P ) 



Dm = {1m 1 - PmH^HPm) 1 ' 2 on H 2 (C P ), 
D° M = (iM 1 - S*P M H*HP M S)^ 2 on H 2 (C P ) 



(2.8) 
(2.9) 



For later purposes we note that S*D 2 M S 



D M . Next define 



uj = 



H 2 ^ 



u{D M Q M ) 



Hi 



D m Rm 



and uj\ Kci . q: Dm = 0. 



(2.10) 
(2.11) 



From the relaxed commutant lifting theory we know that uj is a well defined 
partial isometry with initial space J- = ItciDmQm- Furthermore, the forward 
shift operator V on L 2 is the Sz.-Nagy-Schaifer isometric lifting of V-. Then as 
a consequence of [SI Theorem 4.3] and the above analysis we obtain the following 
result. 



Proposition 2.4. Let G € L^ xp , and let M. be a subspace of H 2 satisfying the 
conditions in (12. ip . Then the central optimal solution Fm to the j\A-restricted 
Nehari problem is given by Fm = G — where &m € L 2 xp has finite Ai- 
norm, the co-analytic part of * s equal to G_ , and the analytic part <&m,+ 
of $jvi is given by 

$M,+ (A) = Wl (J- Xuj 2 )- 1 D m E. (2.12) 
Here E is defined by (ll.6p . and uj\ and uj 2 are defined by (|2.f Op and (|2.1ip . 
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It is this central optimal solution Fm we shall be working with. From 
Corollary 4.4 in [8] (see also (TOj Theorem 1.1]) we know that T = IuiDmQm = 
T>m implies that the central solution of (|2.5|) is the only optimal solution to the 
M -restricted Nehari problem. The latter fact will play a role in Section 2] 

3 Statement of the main convergence result 

Let G S L™ Kp , and let H be the Hankel operator defined by the co- analytic part 
of G. In our main approximation result we shall assume that the following two 
conditions are satisfied: 

(CI) H has finite rank, 

(C2) none of the maximizing vectors of H belongs SH 2 , and the space spanned 
by the maximizing vectors of H has dimension p. 

Note that (CI) is equivalent to G being the sum of a rational matrix function 
with all its poles in B and a matrix-valued H°° function. 

In the scalar case the second part of (C2) implies the first part. To see this 
let p = q = 1, and assume that the space spanned by the maximizing vectors 
of H is one dimensional. Let Sv be a maximizing vector of H. Since S is an 
isometry and V-H = HS, we have u ^ and 

||ff|||M| = ||Jf||||S„|| = \\HSv\\ = \\V-Hv\\ < \\Hv\\ < \\H\\\\v\\. 

Thus the inequalities are equalities, and v is a maximizing vector of H. As 
the the space spanned by the maximizing vectors of H is assumed to be one 
dimensional, v must be a scalar multiple of Sv, which can only happen when 
v = 0, which contradicts v ^ 0. Thus the first part of (C2) is fulfilled. Next 
observe that for p = q = 1 the statement "the space spanned by the maximizing 
vectors of H has dimension one" is just equivalent to the requirement that \\H\\ 
is a simple singular value of H, which is precisely the condition used in Theorem 
2 of the Helton- Young paper [14]. 

As we shall see in Section[5]the two conditions (CI) and (C2) guarantee that 
the solution to the optimal Nehari problem is unique. 

For our approximation scheme we fix a finite dimensional subspace A4q of 
Hp invariant under S*, and we define recursively 

M k = KexS*@SM k -i, ken. (3.1) 

Since M.q is invariant under S* , the space A4q is invariant under S, and the 
Beurling-Lax theorem tells us that Mq = 0i?|, where 8 £ H^ e and can be 
taken to be inner. Using this representation one checks that Aik = HpQz k QH^ 
for each k £ N. It follows that M C Mi C M 2 C • ■ • and Vfe>o Mk = H%. 
Furthermore, 

S*M k cM k and KerS*cM k , k £ N. (3.2) 

Note that the spaces Mk = H 2 z k qH 2 , k = 1, 2, . . ., appearing in [14] satisfy 
(JUT} with M = H 2 QqH 2 . 
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Theorem 3.1. Let G G L^ xp . Assume that conditions (CI) and (C2) are 
satisfied, and let the sequence of subspaces {Mk}ke?t be defined by (|3.ip with 
Mo a finite dimensional S* -invariant subspace of H 2 . Let F be the unique 
optimal solution to the Nehari problem for G, and for each k G N+ let Fk be the 
central optimal solution to the Mk -restricted Nehari problem. Then G — F is 
a rational function in H^ p , and for k G N+ sufficiently large, the same holds 
true for G — Fk. Furthermore, \\Fk — F\\oo ~* f or k — > oo. More precisely, if 
all the poles of G inside D are within the disk D r = {A | |A| < r}, for r < 1, then 
there exists a number L > such that \\Fk — -F||oc < Lr k for k large enough. 

Improving the rate of convergence is one of the main issues in |14j , where it 
is shown that for the case when the poles of G inside D are close to the unit 
circle, that is, r close to 1, convergence with Mq — {0} may occur at a slow rate. 
In [T3] it is also shown how to choose (in the scalar case) a scalar polynomial 
q so that the choice Mo = H 2 Q qH 2 increases the rate of convergence. In 
fact, if the roots of q coincide with the poles of G in D r \D , then starting with 
Mo = H 2 QqH 2 the convergence is of order 0(rp ) rather than 0(r k ). In Section 
[6] we shall see that Theorem 13 . 1 1 remains true if r < 1 is larger than the spectral 
radius of the operator V-\ HM ±, and thus again the convergence rate can be 
improved by an appropriate choice of Mo- To give a trivial example: when Mo 
is chosen in such a way that it includes Im H* , all the central optimal solutions 
F/~ in Theorem 13.11 coincide with the unique optimal solution solution F to the 
Nehari problem. 

4 The central optimal solution revisited 

As before G G L^ p and H is the Hankel operator defined by the co- analytic 
part of G. Furthermore, M is a subspace of H 2 satisfying (|2.1[) . In this section 
we assume that ||JTPmS|| < = H-^-PmII' I r other words, we assume that 
the defect operator D° M defined by (|2.9p is invertible. This additional condition 
allows us to simplify the formula for the central optimal solution to the M- 
restricted Nehari problem presented in Proposition 12.41 We shall prove the 
following theorem. 

Theorem 4.1. Let G G L^ p , and let M be a subspace of H 2 satisfying (|2.1|) . 
Assume the defect operator D° M defined by (|2.9I) is invertible, and put 

Am = D° M 2 S*D 2 M . (4.1) 

Then r S p ec (A_A/i) < 1, and the central optimal solution Fm to the M -restricted 
Nehari problem is given by Fm — G — where G L 2 xp has finite M- 
norm, the co-analytic part of <£>m * s equal to G- , and the analytic part of <$>m 
is given by 

^ M . + {\)=IiH{L-\K M )- 1 K M E = N M {\)M M {\)- 1 (Aei), (4.2) 
where 

N M {\) =nff(J -XS^KmE, M m {X)=I-XE*(I-XS*)- 1 A m E. (4.3) 
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In particular, M(A) is invertible for each A £ D. 

The formulas in the above theorem for the central optimal solution are in- 
spired by the formulas for the central suboptimal solution in Sections IV. 3 and 
IV.4 of 0. 

We first prove two lemmas. In what follows Pm and Rm are the orthogonal 
projections of H£ onto M and S*M, respectively, and Qm — SRm- 

Lemma 4.2. Let Ad be a subspace of Hp satisfying (|2.1|) . Then 

Rm = S*PmS, RmS* = S*Pm, Qm = PmS- (4-4) 

Proof. Note that 

{S*P M S) 2 = S*P M SS*P M S = S*P M S - S*P M (I - SS*)P M S. 

Since I — SS* is the orthogonal projection onto Ker S*, the second part of (|2.ip 
implies that P M (I - SS*) = I - SS*. Thus (S*P M S) 2 = S*P M S, and hence 
S*PmS is an orthogonal projection. The range of this orthogonal projection is 
S*A4, and therefore the first identity in (|4.4I) is proved. 

Using this first identity and Pm (I — SS* ) = I — SS* we see that 

RmS* = S* PmS S* = S*Pm ~ S*Pm{I — SS*) — S*Pm- 

Thus the second identity in (|4.4j) also holds. Finally, 

Qm = SRm = (RmS*)* — (S*Pm)* = PmS- 

Thus P~4")) is proved. □ 

Lemma 4.3. Let G 6 L?^. p) and let M. be a subspace of Hi* satisfying (|2.ip . 
Assume the defect operator D° M defined by (|2.9I) is invertible. Then the range 
J- of the operator DmQm is closed and the orthogonal projection of H 2 onto 
T is given by 

Pt = D m QmD m 2 Q* m Dm- (4.5) 
Proof. We begin with two identities: 

DmPm — PmDm, D° m R m = RmD° m - (4-6) 

Since Pm is an orthogonal projection, the first equality in (|4.6p follows directly 
from the definition of Dm in (|2-8|) . To prove the second, we use the second 
identity in (|4.4I) . Taking adjoints and using the fact that Rm and Pm are 
orthogonal projections, we see that PmS — SRm- It follows that D° M is also 
given by 

D° M = ilM 1 - RmS*H*HSR m ) 1/2 - (4.7) 
From this formula for D° M the second identity in (|4.6[) is clear. 
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Now assume that D° M is invertible, and let P be the operator defined by 
the right hand side of (j4.5|) . Clearly, P is selfadjoint. Let us prove that P is a 
projection. Using the second equality in (|4.6I) we have 

P 2 = D m QmD m 2 Q* m D 2 m Q m D° m 2 Q* m Dm 

— D M Q M D° M 2 (R M S* D 2 m SRm)D° m 2 Q* m Dm 

— D m QmRmD° m 2 (S*D 2 m S)D° m 2 R m Q* m Dm- 

Observe that Q M R M = SR 2 M = SR M = Qm- Since D M = S*D 2 M S, it follows 
that 

P 2 = D M Q M D° M 2 Q* M D M = P. 

Thus P is an orthogonal projection. This implies that P>mQm has a closed 
range, and = P. □ 

Proof of Theorem 14.11 Our starting point is formula (I2.12[) . Recall that wi 
and UJ2 are zero on Ker Q* m Dm- From Lemma 14.31 we know that DmQm has 
a closed range. It follows that lo\ — lo\Pj^ and W2 = where Pjr is the 

orthogonal projection of H 2 onto T — IuiDmQm- Using the formula for Pjr 
given by (|4.5I) . the second intertwining relation in (|4.6[) . the identities in (|4.4I) 
and the definition of u in (j2.10j) . (|2.10p we compute 

u\P>M = lo\PjtDm = luiD mQ mD°m' 'Qm-Dm 

= HH RmD m 2 ' RmS* Dm = D° M 2 RmS* D 2 M 
= HH D° M 2 S* RmD 2 m = HHAmPmi 



and 



— W2P?D M — u 2 D M Q M D M 2 Q M D' M 
= D M R M D° M 2 Q M D 2 M — D M A M P M . 



Furthermore, using the intcrtwing relations in (J4T6J) and the second identity in 
(|4.4I) we see that Rm^m = ^mPm- In particular, Am leaves M invariant. 
Let us now prove that r spC c (A^vf ) < 1- Note that 

r S poc(w 2 ) = r spcc (w 2 Pj) = r spcc (D M RMD° M 2 S* D M ) 

= r spec (RMD M 2 S* D 2 M ) = t svcc (RmKm) — Tsvcc^AmPm)- 

Thus r spec (AMPM) < 1> because u>2 is contractive. Since Am leaves M. invari- 
ant, we see that relative to the orthogonal decomposition H 2 = M © M 1 - the 
operator Am decomposes as 



Aa4 = 



Pm^-mPm * 

(I - P M )A M (I - P M ) 



(4.8) 
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Note that (I — Pm)(7 ~ Rm) = (7 — Pm). Using the latter identity, the formulas 
(|2.8p and (|4.7p . and the intertwining relations in (|4.6p . we obtain 

(7 - P m )*m(I - Pm) = (I- PmW - Rm)^m{I - P M ) 
= (I - P M )(I - R M )S*(I - P M ) 
= (I-P M )S*(I-P M ). 

Thus (7— Pm)^-m{I—Pm) is a contraction. Hence r spoc ((7— Pa/()A.m(7— Pm) < 
1. But then (|4.8|l shows that r spec (A J vi) < 1. 

Next, using that A M M C S*M C .M and ImP = KcrS* C M, we obtain 
for each A 6 D that 

$M.+ (A) = Wi(J- Xuji^DmE = uj x D m {1 - \A M P M y 1 E 

= UHA M P M (I - XKmPm^E = UH(I - \K M P M y l k M P M E 
= UH{I - XAm^AmE, 

which gives formula (|4.2I) . 

Finally, to see that (|4.3I) holds, note that AmS = I. Hence Am is a left 
inverse of S. Since E is an isometry with ImP = Ker5*, we have Am — 
S* + A M EE*. Therefore, for each AeD, 

$A4,+ (A) = HJ7(7 - XAm^AmE = Hf7(7 - XS* - XA M EE*y 1 A M E 
= UH(I - X(I - XS*y 1 A M EE*y 1 (I ~ XS*y 1 A M E 

= uh(i - xs*)- 1 a m e(i - xe*{i - xs*)- 1 a m e)- 1 

= N(X)M(Xy 1 . 

In particular, M(X) is invertible. □ 

Remark. From RmAm = ^mPm we see that Am leaves M. invariant. Thus, 
if M. in Theorem 14.11 is finite dimensional, then $m,+ m (I4.2[) is a rational 
function in H pxql and hence ^>m.+ is a rational pxq matrix function which has 
no pole in the closed unit disk. 

Next we present a criterion in terms of maximizing vectors under which 
Theorem 14.11 applies. 

Proposition 4.4. Assume rankPP^ is finite. Then Dm is invertible if and 
only if none of the maximizing vectors of HPm belongs to SH^ . 

Proof. A vector h G is a maximizing vector of HPm if and only if ^ 
h G P>m- Thus we have to show that invertibility of D° M is equivalent to 

v M nSH 2 p = {0}. 

Assume T> M n SH^ ^ {0}. Thus, using the definition of a maximizing 
vector, there exists Sv with v ^ such that ||7PPm5?;|| = 7x||S'u||. Since S is 
an isometry we see that \\HPmSv\\ = jm\\ v \\- It follows that v is in the kernel 
of D° Ml and hence D° M is not invertible. 
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Conversely, assume that V M n SH 2 = {0}. Note that rank (HPmS) is also 
finite. Hence HPmS has a maximizing vector, say v. We may assume that 
= 1. By our assumption the vector Sv is not a maximizing vector of H Pm- 
Hence 

\\HP M S\\ = \\HP M S\\\\v\\ = \\HP M Sv\\ < ||FPm||||5w|| = j M \\Sv\\ = lM - 

Therefore = ^ m I—S*PmH*HP m S is positive definite, and thus invertible. 
Consequently, D° M is invertible. □ 

For later purposes we mention the following. It is straightforward to prove 
that D° M is invertible if and only if the operator 7^4 1 — HPmSS*PmH* is 
invertible, and in that case we have 

A M P M H* = R M H*Vl(j 2 M I - HP M SS*P M H*)- 1 x 

x{ 1m I - HP M H*), (4.9) 
A M E = -R M H*V1{~( 2 M I - HP M SS*P M H*)- 1 HE. (4.10) 

These formulas can be simplified further using the operators Z and W associated 
to the Hankel operator H which have been introduced at the end of Section [TJ 
see (HHD and (TO). Recall that X = TmH. Since K\ G X = Kcrff*, the 
space X is a reducing subspace for the operators 7^/ — HPj^SS* PmH* and 
7^,1 — HPmH* . Furthermore, 

Ax := [i M I - HP M SS*P M H*)\ X = 1m I x - ZWR M W*Z*, (4.11) 
:= (1m 1 - HP M H*)\x = l 2 M Ix - WP M W*. (4.12) 

Note that Am is invertible if and only if D° M is invertible. Using the above 
operators, (|4.9[) and (|4.10l) can be written as 

A m PmW* = R M W*Z*A M 1 Z M , AmE = -R M W*Z*A M 1 WE. (4.13) 

Corollary 4.5. Let G £ L^ p , and let M. be a subspace of Hi satisfying (|2.ip . 
Assume the operator Am defined by (|4. is invertible. Then the defect opera- 
tor D° M defined by (|2.9p is invertible, and the functions Nm and Mm appearing 
in (|4.3I) are also given by 

N M W = N M!l (\) + N M , 2 (\), (4.14) 

iV M ,i(A) = -UHW*(I - \Z*Y X Z* 'Aj*WE (4.15) 

N M M X ) = UH ( T - X S*)-\I - R M )W*Z*A M l WE. (4.16) 

and 

M M (A) - M M ,i(A) + M M , 2 (A), (4.17) 

Mm,i(A) = / + XE*W*(I - \Z*)~ 1 Z*A M l WE (4.18) 

Mm,2(A) = -\E*(I - XS*)- 1 ^ - R M )W*Z*A M 1 WE. (4.19) 
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Furthermore, if r spcc (Z* A^^m) < 1, then Mm,i(^) * s invertible for \X\ < 1 
and 

MmA*)' 1 =1- XE*W*(I - XZ*A^E M )~ 1 Z*AJ 1 WE, \X\ < 1. (4.20) 

Proof. For operators A and B the invertibility of I + AB is equivalent to the 
invertibility of I + BA. Using this fact it is clear that the invertibility of D° M 
follows form the invertibility of Am- Hence we can apply Theorem l4.il Writing 
Rm as I — (I — Rm) an d using (14.13[) we see that (14. 14|) holds with Nm,2 being 
given by (|4.16[) and with 

Nm,iW = -UH(I - XS*)- 1 W*Z*AJ^WE. (4.21) 

The intertwining relation WS = ZW yields (I - XS*)-^* =W*{I- XZ*)-\ 
Using the latter identity in (|4.21l) yields (|4. 1 5[) . In a similar way one proves the 
identities (j4TT71) - ((4TT^)) . 

To complete the proof assume r spcc (Z* A^S^vi) < 1. Then the inversion 
formula for Mm i(A) follows from the standard inversion formula from [H The- 
orem 2.2.1], where we note that the state operator in the inversion formula 
equals 

Z* - Z*Aj^WEE*W* = Z*A M V ( 1 2 M I - ZWR M W*Z* - WEE*W*) 
= Z*A^( 1 2 M I - W(SR M S* + EE*)W*) 
= Z*A M l ( 1 2 M I - W(SS*P M + EE*P M )W*) 
= Z*A-^{ 1 2 M I ~ WP M W*) = Z*A M 1 Z M , 

as claimed. Here we used the second identity in (|4.4I) . and the fact that PmE = 
E, because ImE = KeiS* C M. □ 

5 The special case where Ai = Hp 

Throughout this section Ai = H 2 , that is, we are dealing with the if^-restricted 
Nehari problem, which is just the usual Nehari problem. Since Ai = H 2 , we 
will surpress the index Ai in our notation, and just write D, D°, T>, T>° , A, etc. 
instead of Dm, D M , T>m, T}° m , Am, etc. In particular, 

J=\\H\\, D = (~f 2 I-H*H) 1/2 , D° = (~f 2 I-S*H*HS) 1/2 . (5.1) 

We shall assume (cf, the first paragraph of Section [3]) that the following two 
conditions are satisfied 

(CI) H has finite rank, 

(C2) none of the maximizing vectors of H belongs SH 2 , and the space spanned 
by the maximizing vectors of H has dimension p. 
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Note that the space spanned by the maximizing vectors of H is equal to Ker D = 
T> , where T> is the closure of the range of D. As H^ — Ker S* © SH 2 , we see 
that 

(C2) ^ Hi = Ker D+SH 2 <^=> H 2 = Ker S*+T>. (5.2) 

Here + means direct sum, not necessarily orthogonal direct sum. 

Let Z and W be the operators defined by (|1.8[) and (|1.9I) . respectively, and 

A = j 2 Ix - ZWW*Z* , 3 = 7 2 /a' — WW*. (5.3) 

We shall prove the following theorem. 

Theorem 5.1. Let G € L 2 xp , and assume that the Hankel operator H asso- 
ciated with the co-analytic part of G satisfies conditions (CI) and (C2). Then 
the operator A defined by the first identity in (|5.3[) is invertible and the Nehari 
problem associated with G has a unique optimal solution F £ H^ p . Moreover, 
this unique solution is given by F = G + — $ +7 where G + is the analytic part of 
G and $+ is the rational qxp matrix-valued H°° function given by 

$+(A) = iV(A)Af(A)" 1 , where 

N(X) = -UHW*{I X - \Z*y 1 Z*A- 1 WE, 
M (A) =7+ \E*W*{Ix - \Z*y 1 Z*A- 1 WE, 

Furthermore, r spec (Z* A _1 S) < 1, and the inverse of M(A) is given by 

MiX)- 1 =1- XE*W*(I X - XZ*A- 1 E)- 1 Z*A- 1 WE. 

Here 3 is the operator defined by the second identity in (|5.3[) . 

The fact that condition (C2) implies uniqueness of the optimal solution 
follows from [2]; cf., Theorem 7.5 (2) in [3]. It will be convenient first to prove 
the following lemma. 

Lemma 5.2. Assume H is compact and (C2) is satisfied. Then the following 
holds. 

(i) The operator D° is invertible, and the range of DS is closed and is equal 
to T> . In particular, the optimal solution to the Nehari problem is unique. 

(ii) The subspace KerD = T> of H 2 is cyclic for S. 

(Hi) The operators u>2 — DD°~ 2 S* D and A = D°~ 2 S*D are well-defined and 
strongly stable. 

Proof. We split the proof into three parts according to the three items. 

Part 1. We prove (i). Since H is compact, the selfadjoint operator D has closed 
range and a finite dimensional null space. Thus D is a Fredholm operator of 
index zero. See |T21 Section XI.l] for the definitions of these notions. Note 
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S is a Fredholm operator of index p. Thus DS is also a Fredholm operator. 
In particular, the range of DS is closed, and hence T := DSH 2 = DSH 2 . 
Moreover, 

ind(L>5) = ind(L>) + ind(S) = -p. 

Here ind denotes the index of a Fredholm operator, and we used the fact f [T2l 
Theorem XL3.2.]) that the index of a product of two Fredholm operators is 
the sum of the indices of the factors. On the other hand, since Ker D n SH 2 
consists of the zero vector only, we see that Ker DS = {0}, and hence, using 
the definition of the index, we have p — codim DSH 2 . But DSH 2 C DH 2 = V 
and, by the third part of (|5.2[) . we have codim £> = p Thus T = T>. The latter 
implies that the central solution of (|2.5[) is the only optimal solution of the 
Nehari problem; see the remark made at the end of Section [2j 

Finally, KerDS = {0} and DS has closed range, yields D° 2 = S*D 2 S is 
invertible. This completes the proof of (i). 

Part 2. We prove (ii). We begin with a remark. From (i) we know that that 
D° is invertible. Thus the operators w 2 = DD°- 2 S*D and A = D°- 2 S*D 2 are 
well defined. Clearly, uj 2 D = DA, and hence uj\D = DA k for k = 0, 1, 2, . . .. It 
follows that 

A k+1 = D°~ 2 S*D 2 A k = D°- 2 S*Dui% , k = 0,1,2,.... (5.4) 

Since uj 2 is a contraction, we conclude that sup fc>0 ||A fe || < oo. 

Our aim is to prove that H 2 — Vfclo S k V^~ . Take h 6 H 2 perpendicular 
to VfcLo S k T >± . The latter is equivalent to S* k h being perpendicular to V 1 - for 
k = 0,1,2,..., that is, S* k h e V for k = 0,1,2,.... Recall that the range of 
D is closed, because H is compact. Thus for each k = 0, 1,2, . . . the vector 
S* k h = Dh k for some h k e V. Thus S* k+1 h = S*Dh k . Since D° is invertible, 
Lemma S31 specified for the case M = H 2 tells us that P := DSD°~ 2 S*D is 
the orthogonal projection of H 2 onto T = T> = ImZ?. Thus for k = 0, 1, 2, . . . 
we have 

S* k h = Dh k = DPh k = D 2 SD°- 2 S*Dh k = D 2 SD°^ 2 S* k+1 h = A*S* k+1 h, 

and by induction h = A* k S* k h. Since lim fe ^ \\S* k+1 h\\ = 0, and sup fc > ||A fc || < 
oo, it follows that \\h\\ = 0. Hence h = 0, and we can conclude that V**lo ^V 1 - = 
H 2 . This proves (ii). 

Part 3. We prove (iii). We already know that W2 and A are well defined. We first 
prove that u>2 is strongly stable, that is, lim^oo wf'u = for any v 6 H 2 . Note 
that uj 2 DS = D. Hence u k DS k = D for k = 0, 1, 2, . . .. Since V s - = KerD, we 
have for any nonnegative integers k, I that u>2 +l DS k V ± = uj^SV 1 - = 0. In other 
words, the kernel of oj k includes X k :— V^ = o S l/ T> ± . Let v £ H 2 . According to 

(ii), we have V^L S " V± = E p- Thus p y* v ~> - with yk=H 2 QX k , and since 
W2 is contractive, we find that 

W^vW = \\uJ k Py k v\\ < 11^11^0. 
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Thus 0J2 is strongly stable, as claimed, and the fact that A is strongly stable 



Proof of Theorem 15. li From Lemma 15.21 (i) we know that D° is invertible, 
and the optimal solution is unique. Since the invertibility of D° implies the 
invertibility of A, we can apply Theorem 14.11 and Corollary 14.51 with Ai = Hp 
to get the desired formula for 4>+. Note that Rh^ = 1% and hence in this case 
the functions appearing in (I4.16|) and (|4.19p are identically zero. 

Put T = Z*AE. Next we show that r spcc (T) < 1. By specifying the first 
identity in flHT3) we see that AW* = W*T, and thus A k W* = W*T k for 
each k € N. Since A is strongly stable (by Lemma 15.21 (hi)), we arrive at 
limfe^oo W*T k x = 0. The fact that H has finite rank, implies that the range 
of H is closed, and hence W is surjective. But then (PFVF*) -1 !^ is a left 
inverse of W* , and T k x = (WW*)~ 1 WT k x -)• if k ->• 00. Thus T is strongly 
stable. Since the underlying space X is finite dimensional, we conclude that 



r spec (Z*AE) = r spoc (T) < 1. 

Finally, since r spGC (Z*AS) = r spcc (T) < 1, the invertibility of M(A) for ||A| 
and the formula for its inverse follow by specifying the final part of Corollary 



6 Convergence of central optimal solutions 

Throughout G £ L^ p and H is the Hankel operator defined by the co- analytic 
part of G. We assume that conditions (CI) and (C2) formulated in the first 
paragraph of Section [3] are satisfied. Furthermore, A4q is a finite dimensional 
^-invariant subspace of H^, and Ma, M.%, ... is a sequence of subspaces 
of Hp defined recursively by (|3.1|) . We set Pk — PM k - From the remarks made 
in the paragraph preceding Theorem 13.11 one sees that 



I-P k = S k (I-P Q )S* k , S*P k = P k - 1 S* P k E = E. (ken). (6.1) 



Here E is the embedding operator defined by (jl.6p . 

In this section we will proof Theorem 13. II In fact we will show that with an 
appropriate choice of the initial space Mq convergence occurs at an ever faster 
rate than stated in Theorem 13.11 We start with a lemma that will be of help 
when proving the increased rate of convergence. 

Lemma 6.1. Let Z and W be the operators defined by (jl.8l) and (jl.9|) . re- 
spectively, and put Xq = WA4q C X. Then Xq is Z -invariant of X — Im W, 
and r spec (Zo) < r spcc (Z). Furthermore, let the operators Zq : Xq — » Xq and 
Wq : — > Xq be defined by Z = Z\x and Wq = Hx W, where Ux is the 
orthogonal projection of X onto Xq . Then 



follows immediately from (15.41) . 



□ 



14.51 for the case when M. = Hp. 



□ 



Z k W(I - P ) - U* Xo Z k W (I - P ), k = 0, 1, 2, . . . . 



(6.2) 
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Proof. Since ZW — WS and M.q is invariant under S, we see that X is 
invariant under Z, and thus 7" spcc (Zo) < r spec (Z). From the definition of Zq and 
W we see that ZW Xq = B* Xo Z and ir^ o W (/ - P ) = W(J - P ). Thus 

z fe M/(/ - p ) - z k n Xa w (i - p ) = n^z*w (i - P.), fc = o, i, 2, . . . 

This proves l|6T2"]l . □ 

Assume < r < 1 such that the poles of G inside D are in the open disc 
B r . As mentioned in the introduction, the poles of G inside B coincide with the 
eigenvalues of Z. Thus r spcc (Z) < r. By Lemma \6. 11 r spec (i?o) < r spec (Z) < r. 
In what follows we fix < r$ < 1 such that r spec (Z ) < r Q < r. We will 
show that the convergence of the central optimal solutions F k in Theorem 13.11 
is proportional to r§- 

For simplicity, we will adapt the notation of Section[5l and write 7, A, N and 
M instead of 7^2 , A. H 2 , N H 2 and M H 2 . Futhermore, we use the abbreviated 
notation P k , j k , A k , E k , and A fc for the operators Pvi fc , lM k , ^ Mk , E Mk , and 
^■M k appearing in Section @] for M = M k . 

As a first step towards the proof of our convergence result we prove the 
following lemma. 

Lemma 6.2. Assume conditions (CI) and (C2) are satisfied. Then Afc — > r 2 A. 
and for k £ N large enough A k is invertible, and A^ 1 — s- r 2 A -1 . 

Proof. We begin with a few remarks. Recall that for M in (|2.ip the operator 
Rm is defined to be the orthogonal projection of H 2 onto S*A4; see the first 
paragraph of Section [2j For M = M k we have S*M k = M k ~\ by (|3.1|) . 
and thus M. = M. k implies Rj\4 k = Pfc-i- It follows that the operator Afc is 
given by A k = ^\l x - ZWP k _ 1 W*Z*; c.f., the second part of (PTTT1) . From 
the invertibility of D° Mk we obtain that Afc is invertible as well; see the first 
paragraph of the proof Corollarv l4.5l The identities in (14. 13)) for A4 = M. k now 
take the form 

AfcPfcVT =Pfc_iVrz*A^Sfc, A k E = -P k ^ 1 W*Z*A k ~ 1 WE. (6.3) 

Observe that 7 2 = ||PPfc|| 2 = ||P fc P*|| 2 - r spcc (HP k H*) = \\HP k H*\\. By a 
similar computation j 2 = \\HH*\\. Thus, using (frTTj) and 

| 7 2 - 7fc 2 | = \\\HH*\\ - \\HP k H*\\\ < \\HH* - HP k H*\\ 

= \\HS k (I - P )S* k H*\\ = \\Z k W (I - Po)W^Z* k \\ 
<\\Z k \\ \\H\\ ||(/-P )||||i/*||||^|| = ||H|| 2 ||Z fe || 2 . 

It follows that 7 2 -^ r 2 7 2 . Next, again by (|6.1|) and (|6.2|) . we obtain 
Afc = 7 2 / - ZWP k ^W*Z* 

= ill - ZWW*Z* + ZWS k ^{I - p )s* k - 1 w*z* 
= A + ( 7 2 - 7 2 )/ + PxX^Woil - Po)W Z* k P Xa . 
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Clearly the second and third summand converge to zero proportional to r^ k , 
and thus we may conclude that A k ~^ r 2 A. 

Since A is invertible by Theorem 15.11 The result of the previous paragraph 
implies that for k large enough A k is invertible and lA^H < L for some L > 
independent of k. Consequently A^ 1 — > r 2 A -1 . □ 

Proof of Theorem 13.11 (with Tq -convergence). We split the proof into four 
parts. Throughout k 6 N is assumed to be large enough so that A k is invertible; 
see Lemma \6. 2 1 

Part 1. Let N and M be as in Theorem I5TT1 Put 

N k ,i( x ) = -UHW*(I - \Z*)- 1 Z*A- 1 WE, (6.4) 
M M (A) = 1+ XE*W*(I - XZ*)- 1 Z*A^ 1 WE. (6.5) 

Since the only dependence on k in N kj i and M k ^\ occurs in the form of A k , it 
follows from Lemma |6~21 that 

Mk,i -> r g M and N k ,i -> r 2 N - (6-6) 

Part 2. From Corollary 14.51 we know that 

N k (X) = N k>1 (\) + iV M (A), N k>2 (X) = UHT k (X)Z*A^WE, (6.7) 
M fe (A) = M M (A) + M M (A), M fc , 2 (A) = -XE*T k (X)Z* A^WE. (6.8) 

Here T k (X) = (I - XS*)' 1 ^ - P k -i)W* . In this part we show that M k ,2 -^r 0. 

Using the first identity in (|6.1[) . the intertwining relation ZW = WS, and 
(|6.2|) we see that 

r fe (A) ={I- A5*) _1 5 fc_1 (J - Po)^" 1 ^* 

= (/ - as*)- 1 ^- 1 ^ - p )w^z; k - 1 u Xo . 

Next we use that 

fe-2 

(7 - AS*)" 1 ^- 1 = ^ X j S k ~ 1 ~ j + A fc_1 (J - AS*)" 1 . 
Thus T fe (A) = r M (A) +r fe)2 (A), where 

r fe ,i(A) = (^A^ fc - 1 ^)(/-p )TU *z* fe - 1 n^ , 

j=0 

r fc , 2 (A) = x k ~ x (i - xs*r\i - p )w^z* k - 1 n Xo . 

Now recall that Mo is ^-invariant, and write So = PoSPq = PoS. The fact 
that Mo is finite dimensional implies r spec (,So) < 1. The computation 

(/ - XS*)- 1 ^ - Po)W£ = (I- XS*)-^ -(I- AS , *)- 1 P W * 

= w*(i - xz*)- 1 -(i- xs^Pow;, 
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shows that (/—AS*) 1 (I—Pq)Wq is uniformly bounded on D. Since r spcc (Z ) < 
ro < 1, we conclude that Tk,2 — >r 0. 

Next observe that E*(j2j=o A^S^ 1 ^') = 0, and thus /:*i\i(A) = for 
each k € N. We conclude that 

Af fc)2 (A) = -AS*r fcj2 (A)Z*A- 1 T^S. 

But then — ►ro implies that the same holds true for Mk,2, that is, —*r 
0. Indeed, this follows from the above identity and the fact that the sequence 
A^" 1 is uniformly bounded. 

Part 3. In this part we show that AT^ )2 — > ro 0. To do this we first observe that 

u,HS k ^ j = nv* -1-5 'H = uv^^PxW = up x z k - 1 - j w. 

Post-multiplying this identity with I — Pq and using (|6.2[) yields 

n/fs^ 1 -^/ - p ) = up Xo z*- l - j w {i - p ). 

It follows that 

fe-2 

fc-i 



Jv fc)2 (A) = ( £ \mp Xo z k ~ x ~ j )w (i - p )w^z* 

nHr k , 2 {\)Z*A^WE. (6.9) 



From the previous part of the proof we know that —^r 0, and by Lemma 
16.21 the sequence A^~ 1 is uniformly bounded. It follows that the second term in 
the right hand side of (|6.9p converges to zero with a rate proportional to r k . 
Note that for A e B we have 

fc-2 fc-2 oo j 

\\Y,xnp Xo z k -^\\ < Y: \\z \\ k -'- j < £ < t±7' 

J=0 j=0 3=1 r ° 

Since r spC c(^o) < r o < 1, we also have — > ro 0. It follows that the first 
term in the right hand side of (|6.9p converges to zero with a rate proportional 
to rj- We conclude that Nk,2 —>r 0- 

Part 4- To complete the proof, it remains to show that Mr (A) — > To M _1 (A) 
uniformly on B. By similar computations as in the proof of Lemma 16.21 it 
follows that S fc — > r 2 S. Hence Z*A^ 1 E k — > r 2 Z*A _1 S. By Theorem 15.11 we 
have r spcc (Z*A _1 S) < 1. Thus for k large enough also r spcc (Z* A^ Hfc) < 1, 
and Mk,i(X) is invertible on B. From the fact that Mk,i -^ r % M, we see that 
MjT{ — > r 2 A/ -1 , with M^l and M~ x indicating here the functions on B with 
values Mfc^A) -1 and A/(A) _1 for each A G B. In particular, the functions 
are uniformly bounded on B by a constant independent of k, which implies 

/ + M~\M K 2 ^ ro /, (/ + M^Mk^)- 1 ^ ro I. 
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As a consequence 

MZ 1 = (M M + = (/ + M^Mk^-Hd^l ^ ro J • M- 1 = M-\ 

which completes the proof. □ 
Concluding remarks 

Note that the functions Mk } \ and Nk,i given by (|6.4|) and (|6.5[) converge with a 
rate proportional to Tq rather than cf., (16.61) . Consequently the same holds 
true for MZ\- Thus a much faster convergence may be achieved when Nk,\MZ{ 
are used instead of N^MZ . However, for the inverse of M^ t i to exist on D we 
need k to be large enough to guarantee r spec (Z* Ak'E.k) < 1, and it is at present 
not clear how large k should be. 

For the scalar case condition (C2) is rather natural. Indeed (sec the second 
paragraph of Section O for the scalar case condition (C2) is equivalent to the 
requirement that the largest singular value of the Hankel operator is simple. The 
latter condition also appears in model reduction problems. In the matrix-valued 
case (C2) seems rather special. We expect that a version of Theorem 13.11 can 
be proved by only using the first part of (C2), that is, by assuming that none of 
the maximizing vectors of the Hankel operator belongs to SH^; cf., Proposition 
14.41 However, note that in that case the optimal solution of the Nehari problem 
may not be unique. 

Computational examples show that it may happen that the approximations 
of the optimal solution to the Nehari problem considered in this paper oscillate 
to the optimal solution when the initial space Mq = {0}. Although the rate of 
convergence can be improved considerably by choosing a different initial space 
Mq, the same examples show that the approximations still oscillate in much the 
same way as before to the optimal solution. This suggests that approximating 
the optimal solution may not be practical in some problems. In this case, one 
may have to adjust these approximating optimal solutions. We plan to return 
to this phenomenon in a later paper. 

Acknowledgement. The authors thank Joe Ball for mentioning the Helton- 
Young paper [14] to the second author. 
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